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Abstract 

We use a recent approach to soft-gluon resummation based on effective field 
theory to implement soft-gluon resummation at NNNLL order for slepton-pair pro- 
duction in SUSY extensions of the Standard Model. This approach resums large 
logarithmic corrections arising from the dynamical enhancement of the partonic 
threshold region caused by steeply falling parton luminosities. We evaluate the 
resummed invariant-mass distribution and total cross section for slepton-pair pro- 
duction at the Tevatron and LHC, matching our results onto NLO fixed-order 
calculations. As a byproduct, we also study virtual SUSY effects in the context of 
soft-gluon resummation for the Drell-Yan process. 
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1 Introduction 



With the advent of the Large Hadron Colhder (LHC), the experimental investigation 
of TeV-scale physics is now fully accessible. The stabilization of the electroweak scale 
requires new particles and interactions in the TeV range, and supersymmetry (SUSY) 
is one of the most compelling scenarios which achieves such a stabilization. For every 
Standard Model (SM) particle it introduces a corresponding superpartner, whose spin 
differs by 1/2 from that of the SM particle. The non-observation of these particles so far 
requires a (soft) breaking of SUSY, and the superpartners are expected to have masses in 
the TeV range. Therefore, it should be possible to observe these particles through their 
direct production at hadron colliders, and indirectly through their virtual contributions 
to SM-like processes. 

In this paper we analyze the production of a scalar-lepton pair at hadron colliders 
1-0], which we consider in parallel with the classic Drell-Yan production of a lepton pair 
5| in SUSY. These two processes can be considered together because they arise from 
the same hard-scattering interaction: the annihilation of a quark-antiquark pair into a 
virtual photon or Z boson, which then decays into a slepton or a lepton pair, respectively. 
The effect of strongly-interacting SUSY particles enters in the hard-scattering interaction 
of both processes only at the one-loop level through the virtual exchange of squarks and 
gluinos. Both processes are interesting and play an important role: Drell-Yan production 
can be considered as a prototype for other collider processes and, among other things, 
its cross section as a function of the invariant mass of the lepton pair can be used to 
search for new heavy resonances. Sleptons are expected to be among the lightest SUSY 
particles, which means that in many scenarios they decay directly into the corresponding 
SM partners and the stable lightest SUSY particle, giving rise to simple signatures such 
as a pair of energetic leptons plus missing energy. 

Both processes have been studied extensively in the past. The calculation of the 
cross section and rapidity distribution at the next-to-leading order (NLO) in for the 
Drell-Yan process in the SM has been accomplished long ago [sl, while the correspond- 
in g re sults at the next-to-next-to-leading order (NNLO) were obtained more recently 
7H 12l|. A study of SUSY QCD and electroweak corrections at NLO was performed in 
13|. Results for the total cross section for slepton-pair production at NLO in were 



obtained in |14|-|l6l| . The main uncertainties in the theoretical predictions arise from the 
imperfect knowledge of the parton distribution functions (PDFs) and from the trunca- 
tion of the perturbative expansion, which introduces a dependence on the unphysical 
factorization and renormalization scales. The two sources of errors are of comparable 
size. In particular, the uncertainty due to scale variations is smaller than in similar 
production processes involving colored particles in the final states, like top-quark pair 
production. This is because, at the partonic tree level, (s)lepton-pair production is a 
purely electroweak process, and therefore at leading order the uncertainty arise only 
from the variation of the factorization scale of the PDFs. The uncertainty from the 
renormalization scale starts at order and is therefore suppressed. A reduction of the 
scale uncertainties is nevertheless desirable, because having a small error on the cross 
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section and the differential distributions allows one to extract interesting information, 
such as the slepton masses, with better precision. 

Reducing the scale uncertainty by computing higher-order terms in the pertubative 
expansion is a task which soon becomes prohibitive. A good alternative for improving the 
theoretical predictions is to resort to soft-gluon resummation methods, which allow one 



to take into account the dominant contributions of the higher-order terms [17|, ll8| . These 



contributions arise from large Sudakov logarithms, which originate as a left-over from the 
cancellation of virtual and real soft divergences in a kinematical configuration where the 
invariant mass M of the (s)lepton pair is close to the partonic center-of-mass energy, and 
hence there is little energy left for additional real radiation. These logarithms must be 
resummed to all orders to improve the convergence of the perturbative expansion. The 
resummation was first achieved at next-to-leading-logarithmic (NLL) order for the Drell- 



Yan invariant-mass distribution in [17|, ll8|, based on a method involving the solution of 



certain evolution equations in Mellin moment space. It was later extended to the rapidity 



distribution 19|], and to N^LL order in 20|,l2l|. In the case of slepton-pair production, the 



resummation was performed at NLL level for the invariant-mass distribution and total 



cross section 



221 ]. Recently, it was shown that Sudakov logarithms can be resummed 



directly in momentum space using methods of soft-coUinear effective theory (SCET) 



23l | . This approach allows for straightforward operator definitions of the various terms 
which build up the cross section, and it avoids the appearance of Landau poles in the 
resummation procedure. The method was used to perform the threshold resummation 
for the Drell-Yan production up to N^LL order in 2J], and then extended to Higgs 
production and top-quark pair production [251428 . 

The aim of our pa per is to extend previous analyses in various directions. First, we 
extend the results of 2J] for the Drell-Yan production of lepton pairs by including the 
contribution from a virtual Z boson as well as SUSY QCD corrections at order as- In the 
case of slepton-pair production, we extend the results of 2^ by performing the threshold 
resummation at N^LL order. While this has a minor effect on the total correction to 
the differential distributions and cross section, it may be relevant for the theoretical 
uncertainty estimate. Moreover, since our resummation is performed by means of SCET 



methods, which constitute an alternative to the Mellin-moment technique used in [22 
we have an independent check of the results in literature. We present results for the 
differential and total cross sections for the Tevatron and the LHC with a center-of-mass 
energy of 7 and 14 TeV. 

In Section [2] we recall the basic formulas for the differential distributions and define 
the kinematics of the threshold region. We then review in Section [3] the factorization 
of the cross section and the resummation of the threshold logarithms in SCET, pre- 
senting the result specific to the SUSY case. Section H] is devoted to a comprehensive 
phenomenological analysis. We estimate the relevance of SUSY QCD corrections and 
the impact of soft-gulon resummation on the invariant- mass distribution and the total 
cross section for slepton-pair production. We also discuss the uncertainties due to scale 
variations and the errors on the PDFs. Our conclusions are given in Section [51 
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2 Kinematics and factorization at threshold 



We consider the production of a (s)lepton pair with invariant mass M in hadron-hadron 
colhsions, at center-of-mass energy ^/s. The process involves the reaction Ni{Pi) + 
^2(-f2) — ^ 7*/Z°* + X, where X represents an inclusive hadronic final state, followed 
by 7*/Z°* — 7- /"(ps) + ^^^(^4) or l^ips) + ^"^(^4)- We start by focusing on the double- 
differential cross section in the invariant mass = g2 and rapidity Y = ^In of 
the (s)lepton pair in the center-of-mass frame, where q = Ps +Pi- This cross section can 
be calculated up to power corrections in perturbative QCD and expressed in terms of 
convolutions of short-distance partonic cross sections with PDFs: 



(Pa 



^ j dXidX2Cij{Xi,X2,S,M,Hf) fi/Ni{Xi,Hf) fj/N^{x2,IJ.f) , (1) 



dM^dY 

where /i/ is the factorization scale, fi/N{x,^f) gives the probability of finding a parton 
i with longitudinal momentum fraction x inside the hadron N, and 



with f3i = — Amj/M"^ denoting the 3- velocity of the slepton in the / rest frame. 

The hard-scattering kernels Cij are related to the partonic cross sections and can be 
calculated as power series in Ug. At leading order (~ a°) the sum involves only the 
channels [ij) = (qq), (qq), with pi = XiPi, p2 = X2P2- At NLO (~ a^) one has to 
take into account {ij) = (qq), (qq), (qg), (gq), (qg)^ (fl'^)- Here we are interested in the 
evaluation of higher-order radiative corrections near threshold, for which it is useful to 
define the quantities 

T 

r= , z = — = , (3) 

S S X1X2 

where s = X1X2S is the partonic center-of-mass energy squared. The partonic threshold 
region is defined by the limit 2 — >■ 1, in which the dynamics of the process is greatly sim- 
plified. Since the partonic center-of-mass energy is just sufficient to create the (s)lepton 
pair, there is no phase space available for the emission of additional energetic partons. 
The cross section is dominated by the terms which are singular in the 2; — )■ 1 limit, which 
correspond to the virtual corrections and the real emission of soft gluons. Such terms 
only arise for the [qq) and (qq) channels. 

We define the couplings of the (s)fermions to gauge boson i following the notation of 



29[. The relevant chiral currents 

f 

= E (atfl^'h + gi'fn'^''fR) , (4) 
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with Pj;^/^ = I (1 =F 75), involve the couphngs 
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(5) 



"w 
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sin0f — 2so 



Here sq^ = sm9]v, c^^ = cos9w, where 9w is the electroweak mixing angle. We 
consider the possibility of mixing between the third-generation sleptons, introducing the 
mass eigenstates fi, f2 and the corresponding mixing angle 9^. 

The leading contributions to the double-differential cross section arising in the near 
partonic threshold 2; — )■ 1 can be written as 24 



2 ^thresh 



dM^dY 



dz 



^0$^ hf^^ / —C{z,M,mq,mg,^if) 



X 



z 

fq/N^Vre^, fif) fq/N^i^e-^, fif) + fq/N^^e^ , fij) /q/^Valv^e"^, fif) 



(6) 



The coefficients /ig ''^ take into account the photon, Z boson, and 7-Z interference con- 
tributions. For lepton-pair production they read 

^2 n 7.2, , i z2 



1 egig}: + g'li ){gl + ) 1 [gl + gf^ )[gl 

2 l-m|/M2 4 (l-m|/M2)2 



/,Z2.1 

9k 



(7) 



while for the production of a slepton pair of type / 



eqi9l + 9W'' , 1 (9 



q2 
L 



9'r)9'^'"' 



2 fl 



(8) 



The hard-scattering kernel C{z, M, m^, rrig, fif) in contains both SM and SUSY QCD 
corrections. The SM QCD corrections are known to two-loop order j9|. The SUSY 
corrections arising at NLO the latter are given by a vertex diagram with a gluino-squark- 
squark loop plus external-leg corrections. They have been calculated in 13, 15, 22. [iol]. 
We have recomputed these contributions and find agreement with results in the literature. 
Integrating over the rapidity, one obtains from ([6]) the single-differential cross section 



da 



thresh 



(T0 2_^ 

q 



1 

h^^ I -C{z,M,mg,mg,fif)ff{T/z,fif) 



where 



^ dx 



X 



[fq/Ni{x,^^f) fg/N2{y/x,fif) + (g o q)] 



(9) 



(10) 



defines the parton luminosity function. 
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3 Factorization in SCET 



The hard-scattering kernel C{z, M, rriq, rrig, fif) depends on the invariant mass M and on 
the masses of squarks and gluinos, rriq and rrig. We will assume in our analysis that these 
three scales are of similar order. For z 1 one can then distinguish two well separated 
mass scales, the "hard" scale fXh ^ M and the "soft" scale /i^ ~ M{l — z)/^/z = \/I{l—z), 
which corresponds to the energy of the emitted soft gluons. The presence of these two 
scales is reflected in the factorization of the coefficient C(2;, M, m^, m^, /ij) into a hard 
and a soft function, 



C{z,M,mg,mg,fif) = H{M,mg,mg, ^f) S{V§{1 - z),fif) . (11) 

Choosing the factorization scale /i/ in ([6]) close to fih or /i^ unavoidably causes the ap- 
pearance of large logarithms in one of the two factors. Threshold resummation addresses 
the problem of resumming these large logarithms to all orders in perturbation theory. 
Recently, the factorization formula has been obtained in the context of SCET 



2J], which is the approach we follow in this paper. It is not necessary to repeat the 
derivation here. It suffices to remember that, considering for simplicity the Drell-Yan 
lepton-pair production through the exchange of a virtual photon with total momentum 
g^, one starts with the expression of the cross section written in terms of the currents 
defined in (jl]), 

d^=^^-0-) J d'xe-^'^-^NMNM\{-9,^)Ji;J^{x)r^^^^^^^ (12) 

The factorization theorem is derived by matching the product of currents onto operators 
in SCET and decoupling soft and collinear interactions by means of field redefinitions. 
In the first step, the current J^^ is matched onto 

Jgr/ Cy{-q^ - ie, nig, m-g, /i) ^ {glXKESh^PLSnXhc + 9RXl^Sh^'PRSnXh^ , (13) 

where the effective fields Xhc = ^Ic^hc and xKc — ^j^hc gauge-invariant com- 

binations of (anti-)hard-collinear quark fields and Wilson lines in SCET. The matching 
coefficient Cy depends on the time-like, hard momentum transfer = —q^ = — M^. 



It is given by the on-shell massless quark form factor [23|, which in the present case 
must be calculated including SUSY QCD corrections. Given that the cross section for 
(s) lepton-pair production involves the current squared, it follows that 



H{M,mq,m-g,iih) = \Cv{-M - ie,mq,mg, ^h)\ ■ (14) 

The soft function arises in a second matching step, in which the soft fields contained in 
the Wilson lines S'„ and Sn in (fT3|) are integrated out. Since it is insensitive to short- 
distance physics, the soft function is the same for Drell-Yan and slepton-pair production. 
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It is defined in terms of tlie vacuum matrix element of a correlator formed by the product 
of the soft Wilson lines contained in the two currents 24 . 

The resummation of threshold logarithms can be achieved directly in momentum 
space by solving the renormalization-group (RG) equations for the hard and soft func- 
tions. In this way, one obtains for the hard-scattering kernel the compact expression 



C{z,M, mq,mg,fif) 



(1 - ^)l-2r, 



sdy In 



r(2r/) 



(15) 



where the definitions of the evolution function U, the coefficient rj and the Laplace- 
transformed soft function sdy can be found in 2J]. The above expression accomplishes 
the resummation of the leading terms for ^ — >■ 1 to all orders of perturbation theory. 
In the SM, the expression for Cy has recently been derived up to three-loop order 
Including the additional virtual corrections of SUSY particles, the perturbative 
expansion up to order 0{a^) can be written as 



31|,l32 



47r 



-M\mlml) 



where c"y and Cy' include the one- and two-loop QCD corrections present in the SM, 
while CysusY ^"^^ ^fsusy represent the additional QCD corrections arising in its SUSY 
extension. In the following we will neglect the two-loop SUSY contribution, which we 



.(2) 



A-n 



-M\mlml 



f^h) , (16) 



assume to be negligible since already c 



(1) 

y,susY 



will turn out to be very small. The explicit 



expressions for Cy^ and Cy' can be found in [2j], while 



.(2) 
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I 2 2\ 

\m^,mg) 



For simplicity we assume degenerate squark masses for R with q = u,d, s, c, b. The 
loop functions Jb and fc are provided in the Appendix. The evolution equation for Cy 
in (fT6|) is the same as for the corresponding coefficient in the SM. As a result, CysusY 
does not depend on the renormalization scale. Note that the Wilson coefficient Cy is 
the same for all currents in (HI). 



4 Systematic studies and phenomenology 

We now present a detailed numerical analysis of the invariant-mass distribution and total 
cross section for slepton-pair production at the Tevatron and LHC. As a byproduct, we 



6 



will also study the rapidity distribution for the Drell-Yan production of lepton pairs. 
Our goal is to estimate the impact of soft-gluon resummation and the relevance of SUSY 
contributions to these observables. To this end, we will either focus on the physical cross 
sections directly or consider the i^-factor defined as 



da 
dM^ 



K{M\mlmlT) 



da 



dM^ 



(18) 



LO 



The theoretical predictions depend on various input parameters, whose numerical values 
are a^miMz) = 1/128, sinOw = 0.23143, Mz = 91.188 GeV, and Tz = 2.4952 GeV. Our 
assumptions for the masses of SUSY particles will be discussed in more detail below. For 
the systematic analyses in Sections O and O we use the PDF set MSTW2008NNLO 
33l . l34l | and as{Mz) = 0.117 with three-loop running in the MS scheme. Using a fixed 



set of PDFs helps to illustrate more clearly the behavior of the perturbative expansion 
of the hard-scattering kernel in higher orders of perturbation theory. 

For appropriate choices of the hard and soft scales fih and /i^, the expressions for 
the hard-scattering kernel given in (IT^ resums the leading singular contributions to the 
partonic cross sections in the limit where the partonic center-of-mass energy a/J is close 
to the invariant mass M of the (s)lepton pair. In fixed-order perturbation theory, the 
leading terms correspond to plus distributions of the form 



1-z 



In 



(19) 



For Drell-Yan type processes, it is well known that, after the hard-scattering kernels are 
convoluted with the parton luminosities, these terms provide the dominant contributions 
to the perturbative series for the cross sections, typi cally accounting for more than 90% 

1). In realistic cases where 



(see e.g. 



24, 26 



of the one- and two-loop corrections 
the ratio r = M^/s is not very close to 1, the dominance of the region 2; < 1 in the 
calculation of the cross section arises dynamically, due to the strong fall-off of the parton 
luminosities. Below we will perform the resummation of the leading terms at different 
orders in RG-improved perturbation theory. At NNLL order, one evaluates (fT6|) using 
the one-loop approximations for the matching functions Cy and sdy along with two-loop 
(three-loop) expressions for the (cusp) anomalous- dimension functions. At N'^LL order, 
one uses the two-loop approximations for the matching functions along with three-loop 
(four-loop) expressions for the (cusp) anomalous dimensions. Note that the two-loop 
virtual SUSY corrections Cyg^gy (113) and the four-loop cusp anomalous dimension 
are currently unknown, but their numerical impact on the N^LL result is expected to be 
negligibly small. In our N^LL results below, we include the known two-loop corrections 
to the hard and soft functions and the relevant three-loop anomalous dimensions. 

Subleading terms can be added to our resummed expressions by matching them to 
fixed-order perturbation theory. To this end, we define 



da 



N"LL+NLO 



da 



da 



NLO I 



I A'/ 



— da 



N"LL| 



=t^f 



(20) 



7 




fiJM M [GeV] 

Figure 1: Relative contribution of the one- loop correction to the soft function to the cross 
section for slepton-pair production at the LHC {^/s = 7TeV), for different values of the pair 
invariant mass M (left). For each value of M we determine the soft scale by taking the point 
at which the correction is minimal (right). 



The first term on the right-hand side denotes the resummed prediction for the cross 
section, while the second one includes those terms that are sub leading in the z — ?■ 1 limit. 
They are obtained by subtracting the fixed-order expression for the leading singular 
terms, derived by setting all three scales equal in expression f|T5|) . from the complete 
fixed-order result. This difference is then expanded to first order in Ug. Since the 
matching to fixed-order theory is somewhat cumbersome, we will sometimes restrict our 
analysis to the leading singular terms only. This will be mentioned explicitly below. 



4.1 Scale setting 

An appropriate choice of the matching scales and /i^, which enter in the resummation 
formula f|T5|) . is crucial for the reduction of the remaining perturbative uncertainties in 
our calculation. In the spirit of effective field theory, the choice of these scales is driven 
by the requirement that the perturbative expansions of the matching coefficients Cy and 
sdy should be well behaved. Since the SUSY contributions to Cy turn out to be very 
small (see below), these effects play no role in the scale-setting procedure, which therefore 



proceeds in analogy with the discussion for the Drell-Yan cross section presented in [24 
For the hard matching scale, we adopt the default choice /i^'^^ = M. 

The soft scale /i^ is set dynamically by minimizing the effect of the one-loop correc- 
tions to the soft function Sdy under the convolution integral This scale therefore 
depends on the value of r = M"^ / s and on the process under consideration. For the case 
of slepton-pair production at the LHC (with -y/i = 7TeV), the result is shown in the left 
plot of Figure [1] for different choices of the invariant mass M. Our default value for the 
soft scale is chosen such that, for fixed M , the contribution of the one- loop correction to 
the soft function to the cross section is minimized. The value of /i^/M for which this con- 
dition is satisfied is shown in the plot on the right. We observe that for sufficiently large 
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values of M the soft scale is indeed much smaller than the hard scale /i/^ ~ M, indicating 
the relevance of threshold resummation. The corresponding plots for the Tevatron and 
the LHC with ^/s = 14 TeV would look similar. For practical purposes, the values of fig 
as a function of M and s can be parametrized by means of the function 

a'^^ - (21) 

where (a, 6, c) = (1.1,3.6,1.9) for the Tevatron, (1.5,4.8,1.7) for the LHC with ^/s = 
7 TeV, and (1.4,3.6,2.0) for the LHC with = 14 TeV. 

After the matching scales have been set, our results still exhibit a residual dependence 
on the factorization scale fif, at which the PDFs are renormalized. As illustrated in 
Figure H] for the case of the i^-factor for slepton-pair production at the LHC, we find 
that after soft-gluon resummation this dependence is significantly weaker than in fixed- 
order perturbation theory, already at NLL order. For this analysis only the leading 
singular two-loop corrections are considered at NNLO, and this is denoted by a star 
(NNLO*). It follows that there is very little sensitivity to the choice of the factorization 
scale. Below we take ^f^^ = M as our default value. 

In Figure [3l we study the i^-factor for slepton-pair production at the LHC, showing 
results obtained at different orders of perturbation theory. Only the leading singular 
terms are included in all cases. The widths of the various bands reflect the scale uncer- 
tainties inherent in the calculations. For the fixed-order results, they are obtained by 
setting the renormalization scale Hr equal to the factorization scale /ij, and varying the 
common scale between 0.5 and 2 times its default value. We performed an independent 
variation of the two scales as well, but the corresponding uncertainty obtained by taking 
an envelope of the maximum deviation from the default value does not differ appreciably 
from the result in Figure [31 This is because at leading order the cross section depends 
on /ij only, while the dependence on Hr starts at NLO and is therefore small. For the 
resummed results, the error bands take into account the uncertainties associated with 
the scales fit, fJ's, and /i/. They correspond to the envelope of the predictions obtained by 
varying all three scales simultaneously between 0.5 and 2 times their default values. We 
verified that the bands obtained in this way do not differ in a significant way from those 
obtained by varying the three scales individually, with the other two scales held fixed, 
and adding then the three uncertainties in quadrature. It is evident that the convergence 
of the perturbative expansion and the remaining scale uncertainties are greatly improved 
by means of soft-gluon resummation. After resummation the three bands nicely overlap, 
and the scale uncertainties are reduced significantly with each order. On the contrary, 
the fixed-order results exhibit a slower convergence and larger scale uncertainties. 



4.2 Impact of SUSY matching corrections 

The existence of SUSY particles would affect our analysis in two ways. First, if sleptons 
exist and are kinematically accessible at the Tevatron or LHC, these particles can be pair 
produced, and via a measurement of their cross section and invariant-mass distribution 
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Figure 2: Factorization-scale dependence of the ii'-factor for slepton-pair production at the 
LHC in fixed-order perturbation theory (left) and after soft-gluon resummation (right). The 
NNLO* and N^LL-|-NLO results contain only the leading singular two-loop corrections. 
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Figure 3: i^-factor for slepton-pair production at the LHC, at different orders in fixed-order 
perturbation theory (left) and including the effects of soft-gluon resummation (right). Only 
the leading singular terms are included. 



one can address questions about the slepton masses and couplings. Secondly, strongly- 
interacting SUSY particles (squarks and gluinos) can affect the hard-scattering kernels 
for both lepton-pair and slepton-pair production at 0{as), via one-loop SUSY QCD 
corrections to the hard matching coefficient Cy in f lTB]) . This second effect is obviously 
more subtle than the first one. We will now address if and to what extent it will be 
possible to probe for virtual effects of SUSY particles in high-precision measurements of 
the standard Drell-Yan cross section. 

In Figure HJ we compare the one-loop SUSY QCD contributions to the Wilson co- 
efficient Cy with the corresponding QCD contributions arising in the SM, for two rep- 
resentative choices of the squark and gluino masses. The first, rrig = 600 GeV and 
rrig = 750 GeV (parameter point Pi), represents an average value for the squark and 
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Figure 4: Comparison of one-loop contributions to the hard matching coefficient Cy arising in 
the SM (sohd hues) and in its SUSY extensions (dashed hues). 
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Figure 5: Comparison of the scale uncertainty of the Drell-Yan ii'- factor for the LHC (light 
bands) with the maximum possible deviation due to loop effects of SUSY particles (dark bands). 
The left plot refers to fixed-order perturbation theory, while the right one includes the effects 
of soft-gluon resummation. Only the leading singular terms are included in the calculation. 



gluino masses close to the point SPSla' (35[. The second, given by rriq = 1200 GeV and 



rrig = 500 GeV (parameter point P2), represents an alternative scenario with a lighter 
gluino and heavier squarks, inspired by scenarios like SPS2 and SPSS [s^. Although 

these points are not yet ruled out in 



close to current limits on sparticle masses |37H41 



the context of a general SUSY extension of the SM. We observe that for all parameter 
choices the SM loop corrections are at least a factor of about 3 larger in magnitude than 
the SUSY corrections. The real part of the SUSY contributions reaches a maximum 
close to the squark threshold, for M > 2mg, where the virtual contributions acquire an 
imaginary part. For values of M below threshold the SUSY corrections are yet much 
smaller. As a result, we find that the virtual SUSY effects are smaller than the residual 
scale uncertainties in the calculation of the cross sections. This is shown in Figure O 
where we compare the factorization-scale uncertainty (both at fixed order and after re- 
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Figure 6: Drell-Yan rapidity distribution at fixed M = 1600 GeV at different orders in pertur- 
bation theory. In both plots, the left half refers to fixed-order perturbation theory, while the 
right half includes the effects of soft-gluon resummation. The plot on the left refers to the SM, 
while the one on the right shows the effects of virtual squarks and gluinos by the dashed bands. 



summation) of the SM-only Drell-Yan i^-factor with the maximum deviation arising if 
we include the virtual SUSY contributions and scan the squark and gluino masses inde- 
pendently over the range between 400 and 2200 GeV. For each value of M we choose the 
points in the (m^, rrig) plane which yield the largest up- and downward deviations. In the 
figure we consider Drell-Yan production at the LHC with ^/s = 14 TeV. At lower energy 
or for the Tevatron the impact of virtual SUSY effects is even smaller. We observe that 
in fixed-order perturbation theory the maximum possible deviation due to virtual effects 
of SUSY particles is always much smaller than the residual scale dependence. After 
soft-gluon resummation the size of the effects is comparable; however, in view of the fact 
that additional theoretical uncertainties arise from the variation of the PDFs (not shown 
in the figure), we conclude that also in this case the combined theoretical error is larger 
than the maximum possible SUSY effect. 

On the one hand, these findings justify an approximation at N^LL order where in 
the two-loop corrections to the hard matching coefficient Cy in flTB]) one neglects the 
SUSY contribution Cyg^sY compared with the corresponding SM contribution Cy\ This 
approximation will be adopted in our N^LL predictions below. On the other hand, it 
appears unlikely that it will be possible to probe for SUSY effects via virtual corrections 
to the Drell-Yan cross section. To illustrate this latter point, we show in Figure [6] the 
result for the Drell-Yan rapidity distribution in our SUSY model Pi (with rriq = 600 GeV 
and nig = 750 GeV) at a value of the pair invariant mass for which the SUSY effects 
are close to maximal (cf. Figure H]). We consider the LHC with ^/s = 14 TeV and re- 
strict our analysis to the leading singular terms only. In the left plot we compare the 
results obtained at different orders in fixed-order perturbation theory with the corre- 
sponding results obtained after soft-gluon resummation. These distributions refer to the 
SM without SUSY effects. We observe that resummation improves the convergence of 
the perturbative expansion and leads to somewhat smaller scale variations at higher or- 
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ders. The right plot zooms in on the central region of the rapidity distribution and shows 
once again the results obtained at leading and next-to-leading order with and without 
resummation. The dashed lines indicate the shift of the NLO* and NNLL bands due to 
the presence of SUSY particles. In both cases (with and without resummation) the shift 
amounts to a small enhancement of the cross section, which however is only a fraction 
of the residual scale uncertainty indicated by the widths of the bands. Note also that 
the additional PDF uncertainty is not included in the plots. 

4.3 Invariant-mass distribution for slepton-pair production 

A clear signal of SUSY would come from the direct production and detection of slepton 
pairs. In this case, the very weak dependence on the squark and gluino masses can be seen 
as an advantage, since it would make it possible to use a measurement of the slepton-pair 
production cross section to extract the mass of the slepton produced, or alternatively 
to set a limit on the slepton mass from an upper limit on the production cross section. 
The sensitivity of the cross section ([1]) to the slepton mass mj arises from the prefactor 

in ([2]), and from the fact that the peak of the invariant mass distribution scales with 
mj. After the determination of the matching and factorization scales in Section 14. 
we are now ready to analyze the impact of soft-gluon resummation on the invariant- 
mass distribution for slepton-pair production. For concreteness, we will consider the 
production of a pair of scalar leptons Z^, as this case has the largest cross section. For 
the SUSY masses we take mj^ = 180 GeV, nig = 600 GeV, and niq = 750 GeV. 

Our results for the slepton invariant-mass distributions are shown in Figure [7] for the 
Tevatron (top), the LHC with ^/s = 7TeV (center), and the LHC with ^/s = 14TeV 
(bottom). Contrary to our previous treatment, from now on we consider the PDFs and a.^ 
at the order which is appropriate for the expansion of the corresponding hard-scattering 
kernels. Specifically, we use LO and NLO PDFs for the LO and NLO fixed-order re- 
sults, and NLO and NNLO PDFs for the resummed results at NLL, NNLL-^NLO, and 
N^LL-|-NLO, since in this case the resummed terms include the bulk of the perturba- 
tive corrections appearing at one order higher in ag- Numerically, the NNLL-I-NLO 
(not shown in the figure) and N^LL-|-NLO results turn out to be very close to each 
other, but the scale dependence of the latter ones is further reduced. We thus consider 
the N^LL-|-NLO approximation as our best prediction. The main effect of soft-gluon 
resummation is to increase the cross sections slightly and to improve the convergence 
of the expansion. The resummation effects become more relevant for larger invariant 
masses. For example, at the Tevatron the increase from NLO to N^LL-|-NLO is 7% at 
500 GeV and 13% at 1000 GeV. The corresponding increases at the LHC are around 2% 
for both ^/s = 7TeV and 14TeV, indicating that at the LHC resummation effects are 
less important. 
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Figure 7: Invariant-mass distributions for slepton-pair production at the Tevatron and LHC. 
The plots on the left show fixed-order results at LO and NLO, while those on the right in- 
clude the effects of soft-gluon resummation at NLL and N'^LL-|-NLO. The bands indicate the 
uncertainty associated with scale variations. 
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Table 1: Total cross sections in fb. The first error refers to the perturbative uncertainties 
associated with scale variations, the second to PDF uncertainties. 

4.4 Total cross section 

We obtain the total cross sections by integrating the invariant-mass distributions over M. 
In Table [1] we present results corresponding to the various approximations discussed in 
the previous section. In this section we provide predictions for both the NNLL+NLO and 
N^LL+NLO total cross sections, so that the effect of including higher-order logarithmic 
terms can be seen. In the table the first error refers to the total scale variation, i.e. the 
variation of /i/ for the fixed-order cross sections and the maximum deviation from the 
default value obtained by varying /i/, fih, and fis simultaneously for the resummed and 
matched results. The second error takes into account the uncertainty of the PDFs at 
the 90% confidence level, which is estimated by evaluating the cross sections with the 
40 sets of PDFs provided by MSTW2008. 

The results shown in the first three blocks of Table [T] refer to the production of a 
slepton II with mass nij^ = 180 GeV (SUSY parameter point Pi, with nig = 600 GeV 
and rrig = 750 GeV) at the Tevatron and the LHC. Note the relevance of the NLO 
correction, which amounts to around 40% for the Tevatron and 30% for the LHC. As 
expected, the resummation effects are larger at the Tevatron, where they amount to 
a 7% enhancement of the NNLL-I-NLO cross section compared with the NLO result. 
At the LHC the resummation gives a smaller 3% additional contribution to the total 
cross section. The additional contribution of the N'^LL+NLO result compared to the 
NNLL-I-NLO approximation is small, below 1%, but performing the resummation at 
N^LL order helps to further reduce the scale uncertainty. 

Since the effect of soft-gluon resummation becomes more important for higher in- 
variant masses, in Table [1] we provide also the total cross section for a heavier slepton 
Il with mass mj^ = 360 GeV (SUSY parameter point P2, with niq = 1200 GeV and 
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Figure 8: Total cross sections for slepton-pair production as a function of the slepton masses. 



rrig = 500 GeV). We only show results for the LHC with a/s = 14TeV, because given 
the small cross sections it would not be possible to observe the production of such heavy 
sleptons at the Tevatron or during the low-energy phase of the LHC. 

In Figure [HI we show the matched N^LL+NLO total cross section as a function of the 
slepton mass. We now consider different types of sleptons, Il^r (with / = e,/i) and fi^2- 
For the staus, we assume a mixing angle 6f = 70°. The cross sections fall off steeply with 
the slepton masses. At the Tevatron it will be difficult to observe sleptons with masses 
exceeding about 250 GeV, while at the LHC it should be possible to observe slepton-pair 
production up to masses in the range 300-400 GeV. In the upper plot we show the cross 
sections at the Tevatron, focusing on the low mass region. The lower plots refer to the 
LHC at 7 and 14 TeV of center-of-mass energy. At the same value of slepton mass, the 
slepton II has the larger cross section, while the slepton Ir has the lowest cross section. 
The cross sections for the production of staus lie in between. We observe that owing to 
the small scale uncertainty at N^LL+NLO, it would be straightforward to extract the 
masses of the sleptons from measurements of the corresponding total cross sections. 

We have compared our predictions for the LO and NLO fixed-order total cross sections 



for slepton production with results provided by the program PROSPINO [15[, finding 
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agreement. Despite the fact that it is difficult to compare our resummed resuhs with 



those presented in |22|, because contrary to these authors we do not consider squark 
mixing, we still find a reasonable agreement. We emphasize that the method developed 
here allows us to resum soft gluons up to the N'^LL order, so that we are able to get a 



smaller scale uncertainty compared with (22 



The main result of our analysis is that, using soft-gluon resummation techniques, we 
can reduce the theoretical uncertainty related to scale variations below the percent level, 
making it a subdominant source of error. This is evident from the results collected in 
Table [H which show that the error due to uncertainties in the PDFs becomes dominant 
beyond NLO. At this level of precision, one may ask whether the NNLO subleading terms 
could also become relevant. While the full NNLO corrections have not yet been calculated 
for the case of slepton pair production, we can estimate their relevance by considering the 
case of Drell-Yan production of lepton pairs, for which the NNLO corrections are availble 
jol. In this case, we find that the additional correction due to the NNLO subleading terms 
amounts to at most 1% of the NLO cross section. It is thus of the same order as the 
scale uncertainty we find at N^LL+NLO. 



5 Conclusions 

If SUSY exists, the pair production of sleptons is an interesting discovery process at 
hadron colliders. Given the simple signature, which for a wide range in SUSY parameter 
space would consist of a pair of energetic leptons plus missing energy, it should be 
relatively easy to detect this process despite its cross section being an order of magnitude 
smaller than corresponding QCD processes, such as squark and gluino pair production. 
We have analyzed the slepton-pair production cross sections at the Tevatron and the 
LHC, along with the related cross section for the Drell-Yan production of a lepton pair, 
with the aim of obtaining accurate predictions by taking into account the effects of soft- 
gluon resummation. This was done by using methods of effective field theory, which 
allow us to perform the resummation directly in momentum space. The factorized cross 
sections in the partonic threshold region are expressed in terms of Wilson coefficiens 
of SCET operators. Solving the RG equations obeyed by these operators allows us to 
resum the large logarithms arising due to soft gluon emissions to all orders in the strong 
coupling constant. 

We have extended the results available in literature in various directions. For the 
Drell-Yan process, we have calculated the effect of virtual SUSY QCD corrections at one- 
loop order. In the case of slepton-pair production, we have extended previous results 
by performing the resummation up to the N^LL level. Moreover, given the fact that 
we perform the resummation not in Mellin moment space but directly in momentum 
space, our results constitute an independent estimation of soft-gluon effects. We have 
provided a detailed phenomenological analysis, presenting results valid for the Tevatron 
and the LHC with -y/i = 7 and 14TeV. We find that the SUSY QCD corrections due 
to the exchange of squark and gluinos are very small in comparison with the SM QCD 
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corrections, and also with the uncertainty in the perturbative calculations. It would 
therefore be challenging to observe the effects of virtual SUSY particles in the Drell-Yan 
rapidity and invariant mass distributions. 

We find that soft-gluon resummation has a small effect on the total cross sections 
for slepton-pair production, ranging from 7% at the Tevatron to 3% at the LHC for a 
slepton mass of 180 GeV. This is a consequence of the fact that resummation effects are 
important only for large values of the invariant mass, corresponding to a region where 
the invariant-mass distribution is very small and gives a tiny contribution to the total 
cross section. Resummation is therefore more important for colliders where r = M^/s is 
larger, i.e. the Tevatron, or for higher slepton masses. On the other hand, it still proves 
useful for the reduction of the theoretical uncertainty due to scale variations. We find 
that the scale uncertainty reduced by about a factor of 2 when going from the fixed-order 
NLO calculation to the N^LL-|-NLO result. The dominant uncertainties then arise from 
other sources, such as the imperfect knowledge of the PDFs, the parameters of the SUSY 
spectrum, Monte-Carlo modeling of the experimental acceptances for these SUSY final 
states, etc. Therefore, the improvement given by soft-gluon resummation constitutes 
only a first step towards a better derivation of the slepton masses from measurements of 
the total cross section. 
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Appendix 



Here we provide explicit expressions for the loop functions Jb and fc appearing in f lT7|) . 
distinguishing two kinematical regimes. Below the squark production threshold, i.e. for 
< 4m?, the two functions are real, while above threshold they develop an imaginary 
part. We first give results for the function fs- Denoting x = 4m?/M^, we obtain 

1 



fB{M'^,m^) = — 1 arctan 
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To express the function fc in a compact form, it is convenient to define 
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while for > 4m? we get 
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In the special case of equal masses, rrig = rriq, these results simplify significantly. We 
then obtain 
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where (with x = Am^/M^ as before) 
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